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\Q , Abstract 

Exact analytic solutions for a class of scalar-tensor gravity theories with a hyperbolic 
scalar potential are presented. Using an exact solution we have successfully constructed a 
D I model of inflation that produces the spectral index, the running of the spectral index and the 
amplitude of scalar perturbations within the constraints given by the WMAP 7 years data. 
The model simultaneously describes the Big Bang and inflation connected by a specific time 
delay between them so that these two events are regarded as dependent on each other. In 
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Q . solving the Fridemann equations, we have utilized an essential Weyl symmetry of our theory 
^ I in 3+1 dimensions which is a predicted remaining symmetry of 2T-physics field theory in 
^ ■ 4+2 dimensions. This led to a new method of obtaining analytic solutions in IT field 
theory which could in principle be used to solve more complicated theories with more scalar 
^ ' fields. Some additional distinguishing properties of the solution includes the fact that there 
■ are early periods of time when the slow roll approximation is not valid. Furthermore, the 
infiaton does not decrease monotonically with time, rather it oscillates around the potential 
minimum while settling down, unlike the slow roll approximation. While the model we 
used for illustration purposes is realistic in most respects, it lacks a mechanism for stopping 
inflation. The technique of obtaining analytic solutions opens a new window for studying 
inflation, and other applications, more precisely than using approximations. 

PACS: 98.80.-k, 98.80. Cq, 04.50.-h. 

Keywords: Big Bang, inflation, infiaton. General Relativity, Weyl symmetry, 2T-physics. 



^ Work partially supported by the US Department of Energy, grant number DE-FG03-84ER40168. 



1 



I. INTRODUCTION 



Most inflation theories involve one or more scalar fields which are called inflatons 
llQQlQ. The slow-roll approximation is a standard technique used in the study of in- 
flation generated by different inflaton potentials. However, for the slow-roll approximation 
to be valid, the shape of the inflaton potential has to be shallow. This is because only with 
a shallow inflaton potential the kinetic energy of the inflaton can be neglected compared 
to the potential energy. Therefore the shape of inflaton potentials have been restricted in 
the past to apply the slow-roll approximation. This approximation cannot be used if we 
want to flgure out the dynamics of inflaton flelds in regions where the kinetic energy cannot 
be neglected. To obtain solutions in these regions, one has to solve the full second order 
coupled nonlinear differential equations. In general, these kinds of equations are difficult to 
solve and are approached by numerical methods. 

In this paper, we will analytically solve a scalar-tensor theory with a scalar field a {x) 
minimally coupled to gravity. The full action of our theory is a standard scalar-tensor theory 
of gravity with a single scalar 

S = Jd'xV^I^^R- \g^'d,ad,a -V{a)^ (1) 

For the present application, a will be the infiaton, but in other applications of our method 
it may have other interpretations. The potential is 

V{a) = (^csinh^ (^y^aj +6cosh^ [^')) 

where h and c are dimensionless free parameters of the potential, and ^ is the reduced Planck 
mass m„ = -, where m„ = = 2.43 x lO^^GeV. For our illustrative model we will further 
specialize to c = 646, and fix some integration constant to a specific value {6 = 0, see below), 
for no particular reason other than plotting some graphs. With this choice of parameters 
we will obtain a model of infiation that matches all the current observational constraints 

B30- 

Our main purpose here is to illustrate the method of analytic computation, as well as the 
underlying important ideas inspired by 2T-physics that lead to this method and potential. 
The essential ingredients that permit us to obtain an exact solution is a reformulation of 
the above theory as a gauge fixed version of a theory with two conformally coupled scalars 



0, s. This structure includes an essential local scaling symmetry, or Weyl symmetry, that 
reduces the theory to one degree of freedom while also generating the Newton constant from 
an initial scale invariant theory with no dimensionful constants. This symmetry structure 
in 3+1 dimensions is a direct outcome of 2T-gravity from which the Weyl symmetry 
emerges as a left-over of general coordinate reparametrizaton symmetry in 4+2 dimensions. 
The surviving Weyl symmetry can be gauge fixed in different ways. One gauge choice gives 
the theory defined above which appears very difficult to solve. Another gauge choice, leads 
to fully decoupled equations for two scalars, for the special potential given above, for the 
case of any conformally flat metric, such as the Robertson- Walker case. These structures 
and the method of solution will be explained in more detail in section (IIIip . 

The observational constraints on inflation include the value of the scalar spectral index, 
the tensor to scalar ratio and the amplitude of scalar perturbation. Following the definition 
of the potential-slow- roll approximation given by Liddle and Lyth jsljol, the validity of the 
slow-roll approximation for some potential V (cr) requires two slow-roll parameters ey, tjv to 
be much smaller than 1 
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« 1. (3) 



^2 V{a) 

It turns out that our model in Eqs. (!2l[T|) requires the infiaton to evolve in a period of time 
from the Big Bang to well into the period of infiation when the slow-roll approximation does 
not apply. During this time, there are several interesting behaviors of the exact solutions 
that cannot be captured by the slow-roll approximation, including the connection between 
the Big Bang and infiation, and an oscillatory behavior of the infiaton field as it settles down, 
which could not be discussed before. After some time well into infiation, the kinetic energy 
decreases monotonically, thus the equation of state asymptotically approaches w —> —1 a.t 
late times. 

Exact solutions for several different inflation potentials have been reported 
III 12 1 [isl [1^ [isl 16|, but none of them fall into the category we discuss in this paper. 



There are several interesting properties of our model. 

• First, if we trace back the dynamics of the cosmological scale factor a (r) analytically, 
there is a specific time when it vanishes a{TBB) = thereby defining the Big Bang. 
This point in conformal time r = r^^ is the definition of beginning of physical "time" 
t {tbb) = as defined by a comoving observer. 



• Second, our exact solution simultaneously captures analytically both the Big Bang 
and inflation. In particular, inflation does not happen right after the Big Bang, rather 
it is determined by the model that it happens around 10^ Planck times tj ~ 10^ tpianck 
after the Big Bang. Furthermore, unlike usual practice in inflation models, we do not 
artiflcially insert a boundary value for the inflaton a (r), rather it is predicted to start 
from inflnity at the Big Bang and then smoothly connect to the inflation period. 

• Third, unlike most inflation theories where the inflaton fleld a (r) decreases mono- 
tonically during inflation, after it drops from its inflnite value at the Big Bang, it 
oscillates around the potential minimum during inflation. So, by contrast to widely 
used approximations in past analyses 5{], in our model time cannot be exchanged with 
fleld strength, showing again that there is no substitute for our type of exact analytic 
solutions. 

Among these properties, the flrst one is tentative, since it is not clear we can apply a 
purely classical theory to the very early moment of the universe where quantum effects are 
expected to play an important role. But in any case we think that the interesting physical 
picture of the relationship between the Big Bang and inflation may still remain as a physically 
correct feature. On the opposite end, inflation in our model does not end by itself, rather we 
need to introduce other mechanisms to end inflation around 10^ Planck times after the Big 
Bang. The mechanism of stopping inflation requires additional physical features involving 
reheating that may be possible to incorporate naturally in a more complete physical theory. 
For this reason we leave it to future investigations, while in this paper the primary focus is 
on introducing a technique for obtaining exact solutions for the scalar-tensor theory with 
the potential ([2]) and showing that it provides an attractive description of the Big-Bang and 
inflation. 

The paper is organized as follows: in section ([II]) we review the phenomenological con- 
straints of inflation theories. In section flIII|) we explain in detail a technique for obtaining 
analytic solutions which we adopt from our study of 2T-gravity. This enabled us to obtain 
analytic solutions for the standard gravity action in ([1]) . In section flIV|) we construct a model 
of inflation that matches all the observational constraints. In section (jV]) we summarize our 
conclusions and point out some future directions. 
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II. PHENOMENOLOGICAL CONSTRAINTS OF INFLATION MODELS 



The latest observational constraints on inflation models from WMAP 7 years data can 
be found in ref In this section we will briefly review the zeroth order and first order 
properties of inflation theories. 

First, an inflation theory must be able to generate a period of inflation which is defined 
as > where (t) is the scale factor of the flat Friedmann-Roberson- Walker (FRW) 
metric in the Einstein frame 

ds"^ = -dt^ + a| (t) [dx^ + dy"^ + dz^) 

= a| (r) (-(ir^ + dx^ + dy'^ + dz^) (4) 

In the second line we have pulled out a common factor a^^ (r) = a\ (t (r)) and defined the 
conformal time r, with t = t{T) given by the relation aE (t) dr = dt. Using conformal time, 
the derivative with respect to time t can be rewritten as ^ = therefore, the Hubble 

parameter H = [(^f-) /a^] can be expressed as 

where the overdot denotes derivative with respect to conformal time. Second, an inflation 
theory should produce more than 60 e-folds of expansion in order to solve the horizon and 
the flatness problems, which means In ^ tfg^n j > 60, where a^*^ and a^^*" are the scale 
factors at the end and the beginning of inflation respectively. The above two constraints are 
the zeroth order properties of an inflation theory that comes from a purely classical gravity 
theory. 

To include the fluctuations of the fields that cause the primordial perturbation, one should 
match the amplitude of fluctuations at the horizon crossing to the current observational 
anisotropy of the CMB. This is explained in detail in , [isl , [lol , 20 1, while here we will 



outline the procedure. We consider the small oscillations for the perturbed FRW metric and 
the inflaton 



ds^ = a| (r) {-(1 + 2A)dT^ - 2diBdx'dr + [(1 + 27^) Sij + didjHr] dx'dx^] (6) 
(t = (Te (r) + 6a (7) 
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In this paper we will provide exact analytic solutions for the backgrounds qe (t) and aE (t) . 
For the perturbations we choose a coordinate reparametrization gauge in which 6a = 0, and 
concentrate on the relevant curvature perturbation R. To compute the amplitude of this 
perturbation, one first expands the action in Eq.([T]) up to second order and then puts the 
background evolution on shell. The resulting action for the perturbation is 



S, 



(2) 



After defining a rescaled amplitude v = zR, using the factor 



ZiT = 



H{r) 



(8) 



(9) 



that depends on background evolution, the action for the perturbation v (r, x) becomes 

1 



(2) 



drd X 



(10) 



The Fourier transform of v (r, x) defines the mode function Vk (t) in momentum space 

d^x 



V [T,X] 



(2vr) 



3 [akVk{r)e'^-^ + alvl{T)e '^'^ 



(11) 



It satisfies the Mukhanov-Sasaki equation which is obtained by minimizing S", 



(2) 



Vk + 



Vk = 0. 



(12) 



This equation is to be solved along with a physical boundary condition that corresponds 
to choosing a particular vacuum. Then one can compute the power spectrum Pji (k), the 
spectral index (k) and the running of the spectral index n' (k) as follows 



\vk{n (fc))r 



Us (k) - 1 = 



(r. (k)) 

d\n{k^Pn{k)) 



(k) given by Ue (t*) H (r*) = k, 
dris (k) 

n [k) 



(13) 
(14) 



dink dink 
Notice that all the above quantities are evaluated at horizon crossing {k = QeH) which 
defines the time r^, (k) . This is because the curvature perturbation is freezed out when the 
wave length 1/k stretches outside the horizon. In general, ^ (given below) is a complicated 
function of rthat renders the analytic solution to Eq.f ll2p difficult or impossible. 

The solution of Eq. (fT2|) has commonly been discussed for the cases in which the back- 
ground ttE (t) ,aE (t) can be approximated by the de Sitter background and the slow roll 



approximations. The factor | can then be approximated by ^ for which Eq.f lT2|) is solvable 
analytically. However, in our model, the behavior of | is far from this common fit function 
Hence we suggest a different fit function to approximate | in the relevant time period 
close to r ~ r* (A;) so that this equation is solvable analytically. 

The tensor perturbation hij is discussed in a similar way, starting with 



a| {^—dr^ + {bij + hij (r, x)) dx^dx^^ 



(15) 



The action for this perturbation is 

1 



S{2) = I drd^x a| (t) 



ill - {duKjf 



1-3 



(16) 



The Fourier transform of hij (r, x) is (k) efj (k) hi (r) e plus the hermitian conjugate, 
where efj [k) is the spin-two polarization tensor and s denotes the polarization. The rescaled 
amplitude defined by /i| = \aEh% is then governed by the following action for the mode ^% (r) 



(2) 



— / drd^k 



dE\ 



(17) 



In our model the behavior of — is similar to the standard case and can be approximated by 



a.E 



a/r"^, so the Mukhanov equation for fif, (r) can be solved in the usual way. As in the scalar 
case, we can then compute the power spectrum for the tensor, (k) , at horizon crossing. 
Finally, another observational quantity, the tensor to scalar ratio r is defined as 

Pt (k) , , 

The phenomenologically allowed ranges for r and n and n' are plotted in Sljoj. A pure 
power law spectrum should predict zero running n' ~ 0. Finally, the amplitude of scalar 
perturbations is also observed to be around P-ji (k) ~ 10^^ [7]. 

In this paper we construct a model of inflation that satisfies all of the above observational 
constraints. 



III. SOLVING THE THEORY ANALYTICALLY 



In this section we will solve analytically the equations of motion for the action ([T]) . Assum- 
ing the scalar field a is homogeneous in space and using the flat Friedmann-Roberson- Walker 
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metric 
a 



there are two well known independent Einstein equations and one equation for 



21|. 



K 

y 



1 . 2 

(J +V{a) 



2a| 



aE 



K 
1, 



2 r 



a 



1 .2 



4 + + V (a) 







(19) 
(20) 
(21) 



Where prime represents the derivative with respect to a. The first two equations are the 
fiu = 00 and fiu = 11 components of Einstein Equations and the third equation is the 
equation of motion for a, all expressed in terms of the conformal time r. These three equations 
are coupled second order nonlinear differential equations. We will obtain all the exact 
solutions of these equations as displayed in Eq. fl53l) and Eq. (1541) when the potential is as 
given in Eq.([2]). To obtain this solution we will use a technique which we developed in the 
context of 2T-gravity that results in a theory with Weyl symmetry when reduced to IT 
shadow. For this reason we give a brief outline of how this is inspired from 2T-gravity. 



A. The 2T approach to ordinary gravity 



Here we will not discuss 2T-gravity itself 22(] 23(], which is a theory in d+2 dimensions. We 



will only mention the crucial property of this formulation, namely that it has the right mix 
of gauge symmetries to eliminate all ghosts from the 2T fields (including those extra timelike 
components in vector or tensor fields) and yield shadow fields in two lower dimensions that 
are ghost free fields in physical interacting IT field theories in d dimensions. Dualities relate 
the many possible shadow IT field theories that emerge in the process of gauge fixing. For 
our purposes here we concentrate only on the so called "conformal shadow" . The conformal 
shadow (like other shadows as well) in d dimensions captures holographically all the gauge 



invariant content of the 2T-gravity parent theory in d + 2 dimensions |24l |. 

The action for the conformal shadow of pure 2T-gravity yields ordinary IT General 
Relativity in d dimensions with some constraints imposed on it. In particular it contains a 
dilaton and the full action has the form 

Sgrav = j d'^xV^ !^Zd(p^R (g) + ^g^'d^cpdA - ^ (0)| , (22) 



where (x) , g^u (x) are the d dimensional shadows of their higher dimensional counterparts 
and the potential is unique V (0) = Xcp^^^ up to an undetermined dimensionless constant A. 
The standard curvature term with the Newton constant G is not permitted as a consequence 
of the gauge symmetries of the parent theory. So there are no dimensionful constants in this 
theory. In this shadow, due to the predetermined constant (fixed by the gauge symmetries 
in 2T-gravity), 

there is an emergent local scaling (Weyl) symmetry in d dimensions which is a remnant of 



general coordinate transformations in the extra 1+1 dimensions 23|. Since the coefficient 
of R (g) must be positive, the dilaton must have the wrong sign kinetic energy to satisfy 
the Weyl symmetry, so is a ghost. Using the Weyl gauge symmetry the shadow dilaton 
can be gauge fixed to a constant 0o (thus eliminating the ghost which would also have 
been a Goldstone boson after condensation), yielding precisely Einstein's General Relativity 
with an arbitrary cosmological constant Sgrav = / d'^x^/^{(|)lR (g) — X(j)ff'^^'^~'^^}, where the 
condensate 0q must be interpreted in terms of Newton's constant = {IGnGd}^^ . 

Matter fields can be added to 2T-gravity, including Klein-Gordon type scalars, Dirac 
or Weyl spinors and Yang-Mills type vectors, all in d + 2 dimensions. There are special 
restrictions on each one of these, on the form of their kinetic energies, and the forms of 
permitted interactions among themselves and with the gravitational fields. These restrictions 
emerge from the underlying gauge symmetry. 

Here we are only concerned with the conformal shadow of this theory when it includes 
only gravity coupled to scalar fields. These scalar fields are all the elementary scalars that 
one would introduce in a complete theory (thus including the dilaton, infiaton, Higgs, SUSY 
partners, GUT scalars, etc.). The emerging conformal shadow then has the following form 



in the language of ordinary field theory in d dimensions with one time 



22j(23| 



S = I d^x^g I ' 2^ Z.. M « . . I 

' ' +z,{ct^'-T.^sl)R{g)~<P'"^''-'^f[^] ' 

where f {si/(f)) is an arbitrary function of its arguments Si/(j). This shadow automatically 
has the local Weyl scale symmetry, 

9',. = e^'^^^g,., 0' = e(^-^/^)^(^)0, = e(^-'^/^)^(^).., (25) 
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with an arbitrary gauge parameter A (x) , which can be verified directly in d dimensions. In- 
deed the predicted special value of Zd and the form of the potential V (0, Si) = 02d/(a!-2) j 
that scales like V (tcj), tsi) = t2Q!/(d-2)y ^^^^ crucial properties to realize this local sym- 

metry. We emphasize that a Weyl symmetry was not one of the gauge symmetries of the 
parent 2T-gravity field theory in d + 2 dimensions, rather it emerges in the conformal shadow 
in d dimensions as a remnant from the general coordinate symmetry in the extra dimensions 



231]. Thus the Weyl symmetry is a signature of 2T physics. 

In this coupling of gravity to matter, demanded by 2T-physics, there is an interesting 
model independent physics prediction 22| to be emphasized. Every physical scalar Si in the 
complete field theory must be a conformal scalar that couples to the curvature term just 
like the dilaton. But to avoid being ghosts the must have the opposite relative sign in 
the kinetic term. Then in the conformal shadow the curvature term is predicted to take the 
form (0^ — sf)R (g) with a required relative minus sign as shown. Hence the gravitational 
constant must emerge from the condensates of all the scalars, not only the dilaton's. This 
predicts a physical effect, that the effective gravitational constant G ~ (0^ — ^^isf)'^ is 
not really a constant, rather it must increase after every phase transition of the universe 
as a whole. Since the dominant part of each field is the condensate after the phase tran- 
sition, ignoring the small fiuctuations, the effective gravitational constant is approximately 
a constant in between the phase transitions. Thus the Newton constant we measure today 
cannot be the same as the analogous constant before the various transitions occurred, such 
as inflation, grand uniflcation, SUSY breaking, electroweak symmetry breaking, etc.. Of 
course the earlier phase transitions supply the dominant condensates in the sum. 

There is also the curious possibility that G ~ (0^ — J2i ^1)^^ could turn negative if the 
other scalars dominate over the dilaton in some regions of the universe, or in the history 
of the universe, thus producing antigravity in those parts of spacetime. The effects of this 
idea on cosmology is one of the motivations that led us to investigate the solutions of this 
theory, as we will do in the rest of this paper. Interestingly, we found that the Big Bang is 
related to the vanishing of the gauge invariant (1 — sj/cfP') at which point the effective 
G changes sign. The familiar portion of the universe (with positive G) and its evolution is 
described by positive values of the factor (1 — sj/ (fP) starting with zero at the Big Bang 
and staying in the positive region throughout the history of the universe. The value of the 
gauge invariant quantity (1 — sj/cjy^) oscillates in the positive region before it reaches its 



10 



asymptotic value 1. As a result, the vanishing of this gauge invariant quantity determines 
the Big Bang.^ 



B. The model for the Big Bang and inflation 

To solve the differential equations for inflation (fT9l) , (l20l) , (12T!) , and even harder ones, we 
will start with a theory of the type fl2^ that includes two scalar fields 0, s, with a Weyl local 
scaling symmetry. Our strategy is to use the Weyl symmetry to make some convenient gauge 
choices. In one gauge the theory will reduce to the standard infiaton theory in Eqs. fl^fTl) . 
while in another gauge it will reduce to a completely solvable theory. Since each one is a 
gauge choice, the solvable theory is dual to the infiaton theory'^. The infiaton equations 
f|T9|) , fl20|) , (pTj) are then solved by transforming the solution from the fully solvable version. 

The starting point is then the action 



S = J d'x^g Q^M^9^09.0 - \g^^d,sd,s + ^ (0' - s') R - <p'f ) 



(26) 



where we have used Z4 = j^. Here / can be an arbitrary function of |, but will be later set 
J' (t) ~ (f ) + b in the present application to reproduce the model of Eq.([2]). Here </> 
and s are real scalars while has the wrong sign in the kinetic term. This makes a ghost 
degree of freedom. With the special coefficient ^ in the coupling of the scalars fields to the 
Ricci curvature, this action is invariant under the local Weyl transformation in Eq. (|25l) . Due 
to this local Weyl symmetry, we can eliminate the ghost degree of freedom by gauge fixing. 
So the theory is actually ghost free. This Weyl symmetry will play a crucial roll in solving 
the theory ([1]) . Note that under the Weyl transformations the quantity (1 — is gauge 

invariant, so the sign or the zeros of the effective gravitational coupling (0^ — s^) ~ G^^ are 
the same for all gauge choices. 



^ Meanwhile, the behavior of the gauge dependent quantity G ~ (0^ — sf varies according to the 
gauge choice; for example it can be taken to be a positive constant in the Einstein frame, but only in 
spacetime regions where it is positive. 

^ This duality, which is in the form of a familiar gauge transformation in the present case, is a simple example 
of a rich set of dualities among IT physical systems predicted by 2T-physics, but missed systematically 
in IT-physics. Choosing a Weyl gauge in the present context amounts to choosing how to embed d 
dimensions in d + 2 dimensions, thus creating a perspective of how an observer in d dimensions perceives 
some "shadow" of phenomena that occur in d + 2 dimensions. For a recent summary of 2T-physics that 



includes a description of such phenomena see [24 1 
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Varying the action (126!) with respect to all its degrees of freedom, we derive the equations 
of motion 

R^iu — -^RQij-u = (27) 
V^,4,^_4,v(^).^/'(^) (28) 

. 

6 

where /' denotes the derivative of / with respect to its argument ^. The energy mo- 



y2s = \sR + ct)^r(-\ (29) 



mentum tensor T^y is 



6 

T, 



-d^(pdy(p + d^sdys - I {g^yV^ - V ^,d^) (0^ - 



(30) 



+g^, (^\d<p ■ d<P - ^ds ■ ds - 0V [i 

The above equations are valid for all gauge choices of the Weyl symmetry. 

Now we select gauges. To start with, we will choose the so called Einstein gauge where 
the coefficient in front of the Ricci curvature is set to the usual gravitational constant 

1^(^1-4)^^ (31) 

The subscripts E indicate that these fields, including the metric {gE)^^, , are in the Einstein 
gauge. So in this gauge, 0^; is related to the field s^; as 

</'i. = ±(4 + ;|)'- (32) 

The ± plays no role because (pE appears always quadratically. Now, inserting this into the 
action (12^ , and further using the following field redefinition 

.^ = ^sinh(^), (33) 

we find that the gauge fixed form of the action fl26l) reduces to the infiaton action ([1]) 
and relates the general infiaton potential V (a) to the general function f {s/(f)) . This result 
indicates that the physics of ([1]) is completely equivalent to the physics of f l26|) for any desired 
potential. Hence we can use the Weyl symmetric version (l26ll in other gauges to tackle the 
solution of the infiaton theory. 

Now suppose we wish to discuss the case of a conformally fiat spacetime defined by the 
line element 

ds"^ = (x) {7]f,ydx^dx'') (34) 
12 



where 7]^^ is the flat Minkowski metric while the scale factor a (x) is an arbitrary function 
of spacetime. The Robertson- Walker metric (jl]) that we discuss later is a special case where 
the scale factor is a function of only time, but for now we are considering the more general 
spacetime dependence in a {x) . 

The equations above fl27ti29p simplify for the conformally fiat metric. They now contain 
three fields a (x) , (x) and s (x) . Under the Weyl transformations we can form two invari- 
ants, which we can choose as (here we state the results more generally for any d, but for the 
application we use d = 4) 

= a('^/2-iVands = a('^/^-i)s. (35) 

All the equations above, when specialized to the conformal spacetime, can be written in 
terms of only these two gauge invariants. 

To proceed, we will choose a Weyl gauge where the metric is actually fiat. In this gauge 
we have a fiat {x) = 1, while (p fiat, s fiat are still arbitrary but they are equal to the gauge 
invariants defined above 

aflat (x) = 1 : = (j)flat, S = Sflat- (36) 

This shows that by choosing the "fiat gauge" we can extract all the gauge invariant infor- 
mation in the equations of motion fl27ti29p for a conformally fiat spacetime. 
In the fiat gauge, since i? (77) = and ^/—rj = 1, the action ( l26i) becomes 

S = ld^X (i id,<Pflat? - \ {d.Sfiatf - <Pfiir^f (^)) (37) 

where, the conformally fiat metric g^i, is actually fiat rj^^iy. 

Now we note that in this gauge the two scalar fields (pfiat and Sfiat (equivalently the two 
gauge invariants (p, s) will have completely decoupled dynamics if f {s/cf)) takes the special 
form f [s/(j)) = c (s/0)'^^*-'^^^^ + h. This is when the infiaton potential takes the special form 
in Eq.(l2]) when = 4. In this case the action describes two completely decoupled scalars 
(pfiat and Sflat were (pfiat has a wrong sign for its kinetic term 

s = / A Q - i (a,,,,,)' - H^ar^ - c4ir-') . (38) 

The equations of motion fl28ll29p in this gauge are 

nj, ;,j,('i+2)/(d-2) ,-, 2(i {d+2)/(d-2) ,r,n\ 

04>flat = -J^b^'fl^^'" USfiat = J^Cs)i^t"' '■ (39) 
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where □ = rj^^d^dp. In addition to these equations one should also impose the constraints 
that follow from general coordinate symmetry (which has been gauge fixed) that imply, 
through Eg. (1271) ■ that the stress tensor in Eq. (l5U]) specialized to the flat gauge must also 
vanish 

= 0. (40) 

There remains solving the equations in this gauge. It should be emphasized that the de- 
coupling of ipfiat and Sflat is valid for the general conformally flat metric, that is, for any 
spacetime dependence of a (a;) , and therefore for any a^; (x) in the Einstein gauge in the 
standard form of the theory of Eqs.( EP^ (i.e. not only the Robertson- Walker case). 

In solving these decoupled equations we should not forget that the acceptable solutions 
in the flat gauge must still satisfy the gauge invariant requirement that the sign and the 
zeros of (1 — s^/0^) must be the same in all gauge fixed versions. In particular if we wish 
to relate to the Einstein gauge in which (1 — s\/(l)\) = {2zj_K^(l)\)~^ is always positive (see 
Eq.( l32|) ). then the corresponding acceptable solutions in the flat gauge must also satisfy 
(l — s'fi^^/ (p'fi^^^) > 0. The zero in the flat gauge at Sfiat/4'fiat = ±1 may be attained only 
when in the Einstein gauge se and (pE tend to infinity according to their relation in Eq.f l32p . 
According to the relation to the inflaton a in Eq. (133|) the zero can happen in any gauge only 
when the inflaton a in the Einstein frame blows up. We will see through exact solutions 
that this gauge invariant zero corresponds to the Big Bang! 



C. Relating the Einstein and flat gauges 

Notice that the quantity ^ is scale invariant, so it can be expressed in various Weyl gauges 
as follows 

S_ _ Sflat _ se_ 

and by using Eq.(l33l) this allows us to relate a to s fiat /(t) fiat as follows 



(41) 



Sflat I 2 , 

= Se[Se + 

(pflat 



tanh (^y/2zdK<j) ■ 



(42) 



Thus the inflaton field a in the Einstein gauge is related to the flat gauge fields as follows 



a 



In 



iZdK 



4>flat + Sflat 



lat Sflat 



— )■ In 

d=A 2k 



4^ flat + Sflat 



lat Sflat 



(43) 
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The other field in the Einstein gauge is the scale factor aE (t) which appears in the equa- 
tions fll9tl2ip we wish to solve. It is related to the fiat fields by a Weyl gauge transformation 
involving a local parameter A (x) consistent with Eq. fl25p . Hence we can determine A (x) 
through the equation 



1 



SE = e^^-'^/^)^^^)./,.,, <Pe = e(^-'^/^)^(^V/Zat = [si + (44) 
So, consistent with the Einstein gauge of Eq. flST]) . we can write 

= 2z,k' {<P%^, - s%^,) , for (4,, - s},^,) > 0. (45) 
Now, using Eq. (l45|) we can figure out the scale factor in the Einstein gauge as 



where we used a fiat = 1- So, Eqs.f l43|l46|) provide the duality transformation for relating 
the Einstein-gauge fields {aE (x) , cr (x)) to the fiat-gauge fields {(pfiat {x) , sjiat (x)) . In these 
expressions we may substitute the gauge invariants 0, s instead of 4>fiat, sjiat- 

It should be emphasized that this approach works in any dimension, not just four. By 
solving the simpler decoupled equations ( l39l) for {(p fiat (x) , s fiat (x)) , and then using this 
duality (which depends on d), we will obtain the solutions for the much more complicated 
coupled differential equations for the fields {qe (x) , a (x)) for any conformally fiat metric in 
the Einstein gauge, i.e. any gj^^^ = qe (x) rj^^, with arbitrary dependence. Special examples 
of such spacetime metrics include AdS^, AdSd~i x S^, AdS^-n x S"', any maximally 
symmetric space, any conformally flat space including singular ones, etc. , and of course 
the Robert son- Walker expanding universe that we wish to discuss next. It is worth noticing 
that these special conformally flat spacetimes in d dimensions, which can all be mapped to 
our decoupled system, have a hidden global 5*0(^,2) symmetry, since they are obtained as 
shadows of the completely flat spacetime in c? + 2 dimensions in the context of 2T-physics 



25||26|. 



D. Solution of inflation equations 



For the problem of inflation we need to specialize to the much simpler homogeneous flelds 
that depend only on the conformal time r and take d = 4. The solution for the complicated 
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equations fll9ti2ip for ge (t) and a (r) can now be obtained by solving the decoupled simple 
equations {(pfiat (t) , Sfiat (t)) plus our duality transformation fl43|46l) . Then, equations ( 139|) 
take the form 

- 0//at = -4fe0j;,t, -Sfiat = 4:Cs}i^,. (47) 

where the derivatives are with respect to r. A first integral is obtained in the form 

-^ + b<PU = -E„ % + c4,, = i?. (48) 

Using the 00 component of Eq.f l27|) . or Eq. fj40|) . we find that the constants are related 
Eg = E(f, = E, where E is the energy density of Sfiat and {—E) is the energy density of (pfiat, 
so 

% - b<PU = % + = E. (49) 



This implies that the total energy density in the fiat gauge is zero which makes perfect sense. 
If the total energy density were not zero then we could not have a static fiat metric. 

The nature of the solutions can now be ascertained intuitively because (l49l) looks like 
the problem of a non-relativistic particle moving in a quartic potential (upside or downside, 
depending on the signs of b, c) at a fixed energy E (which can be positive, negative or zero). 

More precisely, equations f H7|) are second order nonlinear differential equations. A first 
integral is already given by Eq. (ll9|l . and this leads to first order non-linear differential 
equations, 4' fiat = i-^/S-E + 2b(j)ji^^ and Sfiat = ±^2E — 2csji^^, that are easily integrated. 
For the case E = the solutions are simple. For E ^ the solutions are expressed in terms 
of the famous Jacobi elliptic function cn {z\^) as seen below. The function cn{z\m) with 
a more general label m is a doubly periodic meromorphic function (similar to the cosine 
function cos{ujz)) where m is the parameter which determines the period of the function 



(for more information see 27| and appendix [X])). All the possible solutions for different 
regions of the parameters are listed below, where we have abbreviated cn (z) = cn {z\^^ and 
defined the parameters 



c 



46 

c 



T = 2\cE\* T, 6 = constant. (50) 



Here C is a dimensionless parameter that depends on the ratio of b and c, r is a dimensionless 
conformal time. The purpose of introducing ( and r is just to simplify the expressions in the 
tables below. However, we will still use the dimensionful r in some of the solutions where c 
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or E or both of them are zero. The constant 5 is a relative shift in dimensionless conformal 
time. The origin of the conformal time can be changed arbitrarily because of translation 
symmetry of the equations under r — r + Tq . Therefore Tq has no physical meaning but the 
relative time 6 which corresponds to initial conditions has physical meaning. 

There are four free parameters in the solutions of these two decoupled second order 
differential equations, which are E^, Es, S and tq. We have just argued that tq is not a physical 
quantity that can be fixed arbitrarily, also Eq. fl27|) requires E^ = Eg = E. This reduces 
the number of undetermined integration parameters associated with boundary conditions to 
only the two parameters E, 6. Together with the two parameters b and c that determine the 
shape of the inflaton potential, we have totally four parameters in the theory which can be 
adjusted to fit phenomenological observations. 

The following three tables elaborate on all possible solutions of Eqs. fj47|) . corresponding 
to all possible solutions of the Einstein equations ( IT9ll2T|) in the Einstein frame 



b 


E 


c 


4>flat 


Sflat 


b>0 


E>0 


c> 




l+cniCr) 


i 

2 


(f)^cn(r + 5) 






c = 




l+cn{Cf) 


1 

2 


V2E (t + E^5^ 








c < 




l—cn{(^T) 


1 

2 


, 1 
— * 




i 

2 






l+cniCr) 




1 c 1 








E = 


c> 




1 




no solution 












c = 




1 




So = const 














c < 




1 






1 












±y'-2cT + K5 






E <0 


c ^ 


\E\\ 


1 




no solution 




\ b \ 










c < 




1 






E\j 1 








\ b \ 








c 1 cn(r+5) 





Table 1 - Solutions for & > 0. 
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b 


E 


c 


4>flat 


Sflat 


6 = 


E>0 


c> 


±V2e(t + E^6 




(f)^cn(r) 






c = 


±V2E (t + E^5^ 










c < 


±V2E {t + E^5 




1 r 1 - 

C 1 S 1 4 1-C7i(t) 2 




E = 


c> 


00 = const 


no solution 






c = 


00 = ccmst 


Sq = const 






c < 


00 = const 


1 

±V/-2CT + K(5 




E <0 


all c 


no solution 





Table 2 - Solutions for 6 = 0. 



6 




c 


4>flat 




6 < 


E > 


c> 


E 
b 






f ^ cn (r + 5) 






c = 


E 
b 






±V2E (r + E^ri) 






c < 


1 

1 b 






1 r n i 
IB 4 l-cn(r+<5) 2 

1 c 1 l+cra(r+<5) 




E < 


all c 


no solution 





Table 3 - Solutions for 6 < 0. 



Since all solutions of (jyfiat and sjiat are known, it is straightforward to calculate a and 
aE from the duality (or gauge) transformation in the previous subsection and the results 
listed in Tables 1,2,3. One can check that Eq.( l46|) and Eq.( l43|) are actually solutions to 
differential equations f lTI?]) , f l2U]) and fl^Tl) by direct substitution provided (pfiat and Sfiat are 
one of the possible solutions in Tables 1,2,3 that satisfies Eq.f l49p . Verifying this requires 
only the second derivatives supplied by Eq. (H7|) and the first derivatives supplied by Eq.( H9|) . 
Of course, these are consistent with the properties of the Jacobi elliptic functions listed in 
the appendix. 

We have used conformal time as the evolution parameter, however, cosmologists usually 
discuss time using physical time t that is measured by a comoving observer. To convert the 
conformal time to the physical time, one needs to perform an integral. 
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t= I dt = J aEir) dr. (51) 

This follows from the definition of conformal time. In the example we are going to discuss, 
the physical time will diverge logarithmically as conformal time approaches a critical finite 
value Too- This is because as conformal time approaches this critical value, our solution 



indicates that the scale factor in the Einstein gauge diverges as a^; (r) 



With all 



the information given in this section, we can proceed to build a phenomenological inflation 
model. 

IV. A MODEL OF INFLATION 

In this section, we will construct a model of inflation that matches all of the phe- 
nomenological constraints. The specific solution we are using is the one in Table 1 with 
6 > 0, c > 0, > 0. The solution for this range of parameter space is 



lat 



1 — cn {(r) 
l + cniCr) 



Sflat 



V2 



cn (r + S) 



(52) 



From Eqs. fH6|) and ( 143|) the scale factor in the Einstein gauge and the inflaton can be written 
as (recall ( 



46 I 3 



and r = 2 Ici^h t) 



1 — cn (^r) 



and 



~2 



In 



1 + cn {(t 
1 + -^cn (r + 5) 



^[cn{T + 6)y 



1— cn(f r) 



%cn (r + 5) 



V2 



l+cn(Cr) 
l-cn(Cr) 



T 

2 



(53) 



(54) 



It should be noted that the energy parameter E determines the scale of a^;, so could be 
fixed if we want to normalize in the conventional way, qe {today) = 1, but since we are 
discussing the early universe we will keep it as shown above. 

In the remainder of the discussion we will further specialize to a specific value of the 
available constants by taking 6 = and C = | (oi' 646 = c). Although there is no reason 
d priori to prefer this model, we will take this as an example to show that our solution 
is compatible with a phenomenological model of inflation. A study of phenomenologically 
consistent more general parameter space, determined by numerical analysis is expected in 
following papers. 
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The scale factor a^; as function of r is plotted in Fig.l 
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r 

Fig.l - The scale factor aE [j) in the Einstein frame, from Eq.( l53|l with ^ = ^and 5 = 0. 

The interesting behavior of this particular solution is that at r = tbb ~ 0.92 the scale factor 
is exactly zero aE [jbb) = 0. This defines the origin for comoving time t {tbb) = at the 
Big Bang^. Unlike pure exponential inflation, where the conformal time of the Big Bang 
is equal to — oo, our model gives a finite conformal time for the Big Bang. Of course the 
finite value tbb ~ 0.92 is not physically significant since r can be translated by an arbitrary 
amount as remarked earlier. A finite conformal time for the Big Bang is a general property 
of our scalar-tensor theory ([1]) . 

After the Big Bang, the scale factor increases monotonically to infinity at Tqo ~ 7.4. 
Converting conformal time to physical time t for a comoving observer using Eq. (ISTj) , the 
corresponding physical time is infinite^. This is because the integral ( ISTl) for t diverges 
logarithmically when r approaches Too- 

^ At the moment of the Big Bang tbb^ it also happens that the gauge invariant variable {l — s^/cfP') vanishes. 
The solution in terms of (1 — s^fiatl 'i^'^fiat) allowed to change sign in the flat gauge, but if we insist that 
the physics must be the physics described in the Einstein gauge in terms of the comoving time interval 
< t < oo, then we must confine the solution to only (1 — s'jiat/'l^'fiat) — ^- interesting to 

note that in a more general gauge there are solutions in which the gauge invariant quantity (l — s^/(/)^) 
can go through zero and change sign, thus making a transition to a region with antigravity where the 
effective Newton constant G ~ (0^ — s^) ^ is negative. Of course, applying a classical theory at the Big 
Bang, where (l — changes sign is incomplete, however, this can be a guideline for a more complete 

quantum theory of the Big Bang, with possibly new physics insights. 

^ The area under this curve up to some point r gives the comoving time t (see Eq.(|51|)). So, it is intuitive 
to define t = to correspond to tbb while t =infinity corresponds to r = Too- The solution (j53p for the 
scale factor a^; as function of r is actually a periodic function, so there is more to the curve than is shown 
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Inflation is defined as the region Jp- (a^;) > 0, or equivalently ^ (ft) ^ Using (!53|) 
with C = the infiation region is r > 77 ~ 2.87. The conversion from the conformal time to 
the physical time t (f) depends on the magnitude of the dimensionless variable b. This value 
of b can be determined by requiring the scalar perturbation amplitude to be of order 10~^ 
Q]. This requires b to be around 10~^^. Using this small value of b we find that the time 
for infiation t = fj corresponds to tj- = 1.6 x 10"^*^ seconds or 3 x 10^ Planck times after 
the Big Bang. Therefore, after fitting observation for the amplitude, the time of infiation 
relative to the Big Bang is predicted in our model. At the time tj the energy density of the 
universe is 7.9 x lO^^GeV which is three orders of magnitude smaller than the Planck scale. 
This guarantees the theory can be applied at the time of the beginning of infiation. 

The analytic expression for the Hubble parameter, H = that indicates the expansion 
rate of the universe is 



H 



which is plotted in Fig. 2. 



-3/2 



-'flat ^flat) 



laWflat — S flats flat 




(55) 



Fig. 2 - The Hubble parameter H (r). 

From Fig. 2 we can see that the Hubble parameter decreases from infinity at the Big Bang 
monotonically to an asymptotic value. The Hubble rate does not change much after infiation 



in Fig.l, but since the range for the comoving time is aheady infinite in one quarter period in the space 
of T, the physics of the model in comoving time < < < 00 is aheady represented by the portion of the 
curve shown in Fig.l. 
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starts at Tj ~ 2.87, thus during inflation we can approximate the conformal time as r 
but not before. 



The fact that ^ (^ofj ^ ^ when r > tj indicates that our theory does not have a 
mechanism for stopping inflation. However, hke most of the inflation theories, one follows 
the dynamics of the inflaton to a certain point and claims that the inflaton decays into the 
Standard Model particles. This is called reheating 28|] 29|]. Similarly, in our theory, we will 



assume after the theory produces enough inflation, another aspect of the complete theory 
will take over. The reheating mechanism is not the focus of this paper but certainly it is an 
important problem to investigate in the future. 

From Eg. (1541) . we see that when the scale factor vanishes aE {tbb) = the inflaton must 
blow up logaritmically. Hence we find that the inflaton field a drops from infinity at the 
Big Bang to zero in a finite conformal time and then keeps oscillating around the potential 
minimum eventually dropping to at r = Tqo- Unlike most inflation theories, this eliminates 
the need for declaring some arbitrary initial value for the inflaton. As shown in Fig. 3, this 
kind of inflation theory was never reported before. Furthermore, most inflation theories 
based on slow-roll simply assume that inflation ends before the inflaton reaches its potential 
minimum, but this assumption is not valid as seen in the plot of our analytic solution. 
Hence, within the slow-roll approximation, neither the connection to the Big Bang nor the 
oscillation around the potential minimum could be discovered. 




Fig. 3 - The inflaton field drops from infinity at the Big Bang and then oscillates around 
the potential minimum, becoming zero at t = Too- 
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FIG. 1: Fig. 4 - Exact curve for - (f) , solid; quadratic approximation, dash. 



To calculate the power spectrum P-ji^ the spectral index and the ratio of tensor to 

scalar perturbation r, we have to solve the Mukhanov-Sasaki equation (1121) with a proper 
boundary condition. The solvability of this equation hinges on the properties of the function 
z{t) / z (r) whose analytic form in our model is given by 

z _ dl jaE/H) 

-z - ia,/H) ' ^^^^ 
where &e/H is computed analytically by using the expressions for (r) , (t) in terms of 
the solutions (p fiat, s fiat given above as 



0- yrjlatSflat Sflat(Pflatj 2 \l/2 

■ '""atVflat ~ ^flati'flat 



^flat ^flat 



f\ (57) 



The explicit behavior of | in our model is plotted in Fig.4. It seen that this is very different 
than the usual approximation a/r^, therefore as explained earlier in the paper following 
Eq.( !T3|) . we cannot use the standard approximation to the Mukhanov-Sasaki equation. 

To solve the Mukhanov-Sasaki equation in the vicinity of the horizon crossing r ~ r* (k) , 
as defined in Eq. (IT^ . it is appropriate to approximate the curve ^ with the dotted curve 
shown in Fig.4. This corresponds to the upside down dotted parabola in the figure 

^^-A;o'-«'(r- w)'. (58) 
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We emphasize that the horizon crossing time t^, (k) is down the tail past the peak r* {k) > 
Tmax- Here the peak of the parabola is at time r = Xmax, while — /cq is the value of ^ at that 
time I (Tmax) = —^0- Both of these quantities {r^ax, — ^o) ^^^^ the parameter are 
determined by our analytic solution for | from Eq. flFTI) . Our analytic solution in Eq. flS^ 
depends on three parameters, namely , ^, 6) , but ^ depends only on 5) , which in turn 
fully determine the three parameters of the fitting parabola. The Mukhanov- Sasaki equation 
then takes the form 

Vk +[{k^ + kl) + a' (r - w)'] Vk = 0. (59) 

With this fitting function one can obtain analytic solutions for (t) in the desired neigh- 
borhood of r < {k) in terms of hypergeometric functions, and choose an appropriate 
boundary condition at the peak of the inverted parabola at Tmax- 

The curves shown in the various figures in this paper assume the values of 5 = and 
^ = 1/2 as an illustrative example. For this choice of parameters we find the fit function is 
given by f ^ -0.23 - 1.2 (r - 5.66)^ . 

The details of the computation for Vkij), and the corresponding PTi,ns,n' ,r outlined 
in Eqs. fll3lll4IIT8l) . can be found in [10]. Below we plot the predicted spectral index Us (f), 
running of the spectral index n' (f), as well as constrain the tensor to scalar ratio r (f) , as 
follows (where we recall f = 2 \cE\^^^ t) 



n 




confonnal time 



Fig. 5 - The spectral index Ug (r). At r = 7.05 the value is 0.97. 

The observational data restricts the value of the spectral index Ug to between 0.92 to 1 
at 60 e-folds before the end of inflation, and the tensor to scalar ratio between and 0.5. 
These ranges are the two-sigma region which corresponds to 95% confidence limit from the 
WMAP 7 years data jsl. Since in our theory, inflation has to be truncated by hand, we 
just need to know if there is a time during inflation such that the values of Ug and r reside 
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within this range. If such a time exists, we can let inflation develop for 60 more e-folds and 
truncate the theory. 

From the figures we see that at r ~ 7.05, or equivalently t = 3 x 10^^'' seconds after the 
Big Bang, the value of Ug is around 0.97. 




conformal time 



Fig. 6 - The running of the spectral index n' (r). At r = 7.05 its value is close to zero. 



The analytic expression for n' is not very transparent. It is plotted in Fig. 6. We can see 
that at r ~ 7.05 the running is around 0. This is within the upper bound of the WMAP 7 
years data estimates of around 0.05. 

Meanwhile, the observational upper bound of the tensor to scalar ratio r is satisfied in our 
model and merely constrains the boundary condition of tensor mode perturbation equation 



at Try 



io|. 



V. CONCLUSION AND FUTURE DIRECTIONS 



In this paper we have established a definite relation between the Big Bang and inflation, 
where neither one could be considered as separate events due to separate physical reasons. 
This is a new point of view for the cosmological evolution of our universe. 

We have also presented both a new method for solving exactly certain scalar-tensor 
theories, as well as a sufficiently realistic new model of inflation. The theory has very few 
parameters, but yet it makes phenomenologically consistent interesting predictions based on 
an analytic solution, and simultaneously reveals several new physical features not discussed 
before. 

The method consists of enlarging the action to a form with local Weyl symmetry, with 
some special potential. Both of these aspects were inspired by 2T-physics. We can then 
solve the theory in a gauge where exact solutions are obtained. 

25 



This method should be regarded as one of the useful technical by products of 2T-gravity, 
while the physics content is also one of the novel predictions of 2T-physics. The new features 
of our solution could have been obtained within IT-physics, but was not, thus demonstrating 
that IT-physics lacks the guidance that 2T-physics supplies systematically in the form of 
hidden dualities and hidden symmetries as explained in 2J] 25| 26 1. 

We could in principle use the same techniques to obtain exact solutions for more general 
inflation models with more than one scalar field, and more complicated potentials. As we 
have seen in this paper, a theory that appears intractable may be presented as a complicated 
IT-physics shadow of a 2T-physics theory, while in another shadow that corresponds to a 
more tractable gauge choice, the theory is handled much more easily and even solvable. As 



illustrated in 



24j | 25| 26| 2T-physics offers this kind of new insights that are not available 



in IT-physics. In this paper we have essentially used this duality idea of 2T-physics for a 
specific case in the context of 2T-gravity, but we expect more general applications of this 
concept in future work. 

The exact solutions obtained with our methods allowed us to analyze the theory more 
precisely than using the slow-roll approximation. This revealed properties of the theory that 
could not be found by slow-roll analyses, such as the predictions for the time delay between 
the Big Bang and inflation or the oscillatory behavior of the infiaton. 

Satisfying the constraints coming from the first order phenomenological parameters, the 
spectral index, the running of the specrtral index and the amplitude of scalar power spec- 
trum, we have constructed a phenomenological model, with the action ([T]) and infiaton 
potential 

6^^ 





V{a) = ( -j 6 ^64sinh^ ^ ^^ -a j + cosh^ | \l -a | | , (60) 

where 6 is a positive dimensionless parameter of order 10^^^. A small parameter like h is 
common in general inflation theories. This makes the coefficient of the potential V of order 
(;^)^ 6 ~ (6 X lO^^GeV)^ , revealing perhaps an interesting scale since it is not too far from 
the grand unification (GUT) scale. 

This particular theory has several interesting behaviors. However, it does not predict 
when inflation ends, which is one of its weak points. In a complete theory, a modification 
around the GUT scale mentioned above could provide the desired mechanism to end inflation. 

It is also interesting to note that we obtain a cyclic cosmology, not unlike ref.jsol, under 
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the following conditions 

1. A special value of the integration constant S given by 

the quarter period of the elliptic function cn {f + 6) . (61) 

2. A quantized value of the parameter C = f ? with integer n. 

Under these conditions our solution describes an universe that expands from a big bang and 
shrinks to a big crunch, and repeats this cycle indefinitely. In the meantime the equation of 

state w 

^ a%/2al - V (g^) 

al/2al + V{aE) ^""'^ 
never goes below —1. In this case the effective gravitational "constant" [(p'jiat (t) — s'ji^^ (r)) 
that appears in the action Eq. ( l26i) never changes sign in any cycle although it becomes zero 
at each big bang or big crunch. This cyclic cosmological scenario will be discussed in more 



detail separately in a future paper [31 1- 



The applications of our equations are not restricted to inflation theories. They might also 
be useful to discuss dark energy to explain the current expansion of the universe. However, 
this will require a fine tuning of our parameters. If one wants to use a similar model for dark 
energy, b has to be as small as 10~^^°. This is because the current Hubble time is way too 
big compared to the Planck time. Including supersymmetry in our approach may play a role 
to suppress the value of emergent cosmological constant at late times (see 32|). A modified 
version of our theory might be extendible up to the current era, and then be applicable also 
to the physics of dark energy. 
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Appendix A: Properties of Jacobi elliptic functions 

The Jacobi elliptic functions used in this paper are sn {z\m) , cn {z\m) and dn{z\m). 
There are two parameters z and m for each of them. The parameter m determines the 



27 



period under translations of the parameter z. The period is given by 

dQ 



•period = 4 * 



r- 



\ ■ 

2 



(Al) 



(l — m sin^ 6*) 

With m = the period for all the Jacobi elliptic functions is 27r. For m = these functions 
reduce to the familiar sine and cosine, sn (-zjO) = sin {z\ cn {z\0) — cos (z), and then dn (z) — 
1. In our solutions we have m = |. This results in the period to be about 7.42 in z, but 

note that z is related by a factor to r or it is translated by d in the various entries in Tables 
1,2,3. The functions sn {z\^) , cn and dn (^||) are plotted in Fig.7 




Fig. 7 Plots of cn(2;|i) (sohd), sn(2;|^) (dash) and dn(2;||) (dash-dot). 

We can see that the behavior of sn and cn are very similar to their counterpart trigonometric 
functions. They also satisfy properties similar to trigonometric functions, such as 



(sn {z\m))'^ + (cn {z\m))'^ — 1; m (sn {z\m))'^ + (dn {z\m))'^ — 1 



(A2) 



The derivative of Jacobi elliptic functions are given in terms of expressions somewhat 
similar to those for trigonometric functions. 

-^sn {z\m) — cn {z\m) x dn {z\m) , 
dz 



d_ 

dz 
d 

dz 



cn {z\m) = —sn {z\m) x dn {z\m) , 
dn {z\m) — —m x sn {z\m.) x cn {z\m) 
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One can look up more properties of Jacobi elliptic functions at [2 
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